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ABSTRACT 

Horizontal  and  vertical  thermal  structures  were  examined 
in  a  region  of  the  Northeastern  Pacific  Ocean  off  Northern 
California.  The  observations  were  acquired  on  three  cruises 
as  part  of  the  ONK-sponsored  Ocean  Prediction  Through 
Observation,  Modeling  and  Analysis  (OPTOMA)  research 
program,  centered  in  a  region  of  the  California  Current 
System,  ca.  37  to  39°N,  124  to  126°W,  during  June  and  July, 
1983. 

The  horizontal  temperature  correlation  scale  was  between 
30  and  50km,  which  was  a  significant  factor  when  comparisons 
were  made  between  measured  horizontal  thermal  structures  and 
those  retrieved  from  Fleet  Numerical  Oceanography  Center* s 
(FNOC)  analyses,  which  had  grid  spacing  of  320km  (TEOTS)  and 
40km  (EOTS)  . 

Operational  (modeled)  analysis  fields  were  not  in  agree- 
ment with  the  observed  fields.  The  major  discrepancies 
occurred  in  the  magnitude  of  the  mixed  layer  depth  and  the 
shape  of  the  horizontal  temperature  fields  (maps)  . 

The  operational  significance  of  differences  between 
modeled  and  observed  thermal  structures  was  assessed  in 
terms  of  their  effect  on  low-frequency;  i.e., less  than  1Khz, 
acoustic  propagation  utilizing  the  FACT9 H  and  PE 
Transmission  Loss  models. 

Mixed  layer  depth  differences  produced  significant  disa- 
greements between  direct  ranges  computed  from  model  and 
observed  temperature  profiles.  The  effect  was  most 
pronounced  at  higher  frequencies  and  when  both  source  and 
receiver  were  shallow;  i.e.,  both  at  20m. 

A  comparison  was  made  between  average  depth/temperature 
profiles  from  July,  1982,  July,   1983,   and  FNOC  climatology 


to  obtain  a  measure  of  the  effect  of  interannual  variability 
in  the  domain.  This  comparison  showed  that  a  significant 
temperature  anomaly  existed  in  the  upper  400m  in  1983 
compared  to  1982  due  to  El  Nino  and  that  this  anomaly  was 
not  represented  by  the  FNOC  climatology. 

The  differences  between  modeled  and  measured  thermal 
structures  are  believed  to  be  related  to  thermal  structure 
model  resolution,  model  sensitivity  to  input  data,  short 
scales  of  spatial  variability  and  non- representative  clima- 
tology for  the  domain. 
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I-  I5X1QD0CTI0N 

A.   BACKGROUND  AND  OBJECTIVES 

The  dependence  of  acoustic  propagation  on  the  ccean 
thermal  structure  and  the  increased  demand  from  the  Navy  for 
accurate  acoustic  forecasts  to  facilitate  planning  and 
economical  use  of  resources  has  led  to  an  increasing 
requirement  for  a  reliable  and  accurate  model  to  represent 
ocean  thermal  structure.  Such  a  requirement  is  not  easily 
met,  however,  because  of  the  limited  number  of  observations 
in  the  typical  ocean  thermal  structure  data  base,  approxi- 
mately 150  XBT  measurements  per  day  in  the  northern  hemi- 
sphere [Clancy  et  al.,1983].  In  view  of  the  upper  ocean's 
spatial  and  temporal  scales  of  variability  such  a  small  data 
base  represents  a  significant  handicap  to  the  performance  of 
any  ocean  analysis  mcdel. 

The  availability  of  an  ocean  thermal  structure  data  base 
from  a  relatively  densely  sampled  area  of  the  Northeastern 
Pacific  Ocean  (approximately  200km  west  of  point  Arena, 
Northern  California)  made  it  possible  to  evaluate  the  accu- 
racy of  the  output  of  two  operational  ocean  thermal  struc- 
ture analysis  systems. 

The  objectives  of  this  research  are  firstly  to  obtain  a 
description  of  the  characteristics  of  the  ocean  thermal 
structure  from  a  pcrtion  of  the  ocean  (ca.37°N  to  39°N; 
124° W  to  126°W)  sampled  during  the  three  cruises  comprising 
CPT0MA5  which  took  place  in  June  and  July  1983.  (OPTOMA  - 
Ocean  Prediction  Through  Observations,  Modeling,  and 
Analysis,  is  a  joint  Harvard/NPS  project  'intended  to 
acquire  field  data  to  characterise  synoptic  scale  eddies 
over   a  domain   in   the  California  Current  off   Northern 
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California,  and  to  ' set-up'  an  eddy-resolving,  statistical/ 
dynamical,  limited  domain,  open  boundary  numerical  ocean 
prediction  model. 'The  '5*  indicates  the  fifth  series  of 
cruises  of  the  project).  Secondly,  the  analysed  ocean 
thermal  structure  is  to  be  used  as  a  'sea  truth*  by  which  to 
evaluate  the  accuracy  of  two  real-time  ocean  thermal  anal- 
ysis systems  presently  in  use  at  the  Fleet  Numerical 
Oceanography  Center  (FNOC)  ,  namely,  the  Expanded  Ocean 
Thermal  Structure  (ECTS)  analysis  and  the  Thermodynamic 
Ocean  Prediction  System  (TOPS) -Coupled  EOTS  (TEOTS)  anal- 
ysis. The  final  objective  is  to  determine  the  operational 
significance  of  deviations  in  the  modeled  thermal  structure 
from  the  •sea  truth'  by  utilizing  the  output  from  opera- 
tional acoustic  models  to  provide  a  quantitative  measure  of 
the  effect  of  these  thermal  structure  differences  on 
acoustic  propagation  in  the  ocean.  Two  low-frequency, 
passive- acoustic  models  were  used  to  provide  a  measure  of 
the  variability  of  acoustic  energy  propagation  which  could 
be  expected  from  the  acoustic  models  alone.  These  then, 
would  provide  a  benchmark  by  which  the  variability  due  to 
differences  between  modeled  and  analysed  thermal  structures 
could  te  judged  significant  or  not.  The  two  acoustic  models 
used  were  the  FACT9H  (Fast  Asymptotic  Ccherrent 
Transmission)  transmission  loss  model  and  the  AESD  (Acoustic 
Enviromental  Support  Detachment)  PE  (Parabolic  Equation) 
model  [ Brock, 1 978  ]. 

B.   STUDI  AREA  AHD  ITS  OCEANOGBAPHY 

The  area  investigated  is  located  in  the  California 
Current  System  (CCS) ,  south  of  the  Mendocino  Escarpment.  The 
region  is  approximately  bounded  by  37<>  to  39°  N  and  124°  to 
1260W.  The  bathymetry.  Fig.  1.1  ,  shows  the  area  to  be 
offshore  of  the  continental  slope,  with  an  average  water 
depth  of  4000m. 


16 


CN2 


O 

O 

OT 

•("H 

d 

o 

«3 

d 

CO 

cd 

*■«    ^* 

Cn  /^ 

ft) 

0) 
M 

-a 


a 
a> 
o 
«j 

■n 
•o 

o 
<d 

a 

•H 

id 

a 
o 

Q 

>i 

S3 
■M 
V) 

<M 
O 

5H 
M 
■M 
0) 

a 

-p 
id 


0) 

en 
•H 


17 


The  CCS  is  composed  of  four  currents:  the  California 
Current,  the  Davidson  Current,  the  California  Undercurrent, 
and  the  Southern  California  Current.  Of  these  four  currents 
the  southward  flowing  California  Current  is  the  major  influ- 
ence in  the  study  domain.  This  current  forms  the  eastern 
boundary  current  of  the  large  anticyclonic  gyre  centred  near 
the  Hawiian  Islands.  As  an  extension  of  the  Westwind  Drift, 
the  California  Current  occurs  between  the  North  Pacific 
atmospheric  high  pressure  system  and  the  semi-permanent 
thermal  low  positioned  over  Central  California.  Contrary  to 
the  classical  picture  of  a  broad,  slow,  shallow  climatolo- 
gical  mean  eastern  boundary  current,  the  CCS  is  known  to 
comprise  mesoscale  eddies,  meandering  jets  and  turbulent 
filaments  [ Mooers  and  Robinson,  1984]  together  with  their 
associated  complex  horizontal  and  vertical  thermal  struc- 
tures. The  mean  current  is  broad  (approx.  1000km  wide), 
shallow  (less  than  500m  deep)  and  sluggish  moving  (approx 
0.9km/hr).  The  southward  to  southeastward  winds,  which 
prevail  during  the  spring  to  the  fall  months  south  of  40°  N, 
produce  upwelling  along  the  coast,  extending  50  to  100km 
offshore.  Upwelling  begins  in  February  or  March  and 
continues  through  August  in  the  study  domain.  The  cold, 
usually  high  salinity  water,  brought  into  the  surface  layers 
by  this  process,  is  then  moved  offshore  by  the  smaller  scale 
perterbations  in  the  CCS,  further  complicating  the  ocean 
thermal  structure.  The  most  intense  upwelling  sites  cccur  to 
the  south  of  capes  and  headlands  [Reid  et  al.,1958]  with  the 
most  conspicuous  centres  of  upwelling  adjacent  to  the  study 
domain  at  35<>  N  and  410  n  [Sverdrup  et  al.,  1942  and  Reid  et 
al.,  1958]. 

Mean  monthly  upwelling  velocities  for  the  west  coast  of 
North  America (Fig  1.2;  from  Heath  (1983)),  based  on  Nelson's 
(1977)  wind  stress  curl  data,  clearly  show  upwelling  to  be 
at  its  most  intense  in  July  off  Northern  California.   Hence, 
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the  maximum  horizontal  and,  because  of  the  differing  water 
masses  involved,  vertical  variations  in  thermal  structure 
should  be  expected  to  occur  in  the  study  domain  during  the 
period  under  investigation. 

In  addition  to  the  seasonal  variations  in  thermal  struc- 
ture, inter-annual  variability  is  also  evident.  Significant 
increases  in  temperature  off  Central  and  Northern  California 
occurred  during  the  winters  of  1972,  1976,  1979  and  1983, 
corresponding  to  El  Nino  episodes  [Breaker,  1983].  Of 
particular  interest  are  the  effects  of  the  major  El  Nino 
warming  of  1982-83.  A  brief  analysis  of  the  possible  opera- 
tional ccnseguences  of  such  an  anomolous  event  is  presented 
in  Chapter  V. 
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Figure  1.2    Bean  Honthly  Opwelling  Indices. 
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II.  DESCRIPTION  OF  COM POT ER  MODELS 

A.   OCEAN  ANALYSIS  MCDELS 

T.   Expanded  Ocean  Thermal  Structure  (EOTS)    Analysis 
i?2^el 

The  Expanded  Ccean  Thermal  Structure  (EOTS)  analysis 
[Mendenhall  et  al.,  1978]  has  been  the  OS  Navy's  operational 
thermal  analysis  system  for  the  past  several  years  and 
provides  ocean  thermal  structure  input  to  most  of  the 
acoustic  predictions  generated  at  FNOC. 

EOTS  is  run  en  a  63x63  grid  for  both  hemispheric  and 
regional  coverage.  The  grid  spacing  for  the  hemispheric 
product  is  200  to  400km,  while  that  of  the  regional  version 
is  typically  40  to  80km.  For  the  SOCAL  region,  in  which  the 
study  domain  resides,  the  spacings  are  320km  and  40km 
respectively. 

The  EOTS  analysis  system  is  based  on  a 
Fields-by-Information-Blending  (FIB)  methodology  [  Holl  and 
Mendenhall,  1971;  Hell,  Cuming  and  Mendenhall,  1979]  and  is 
a  numerical  blending  technique  whereby  the  difference 
tetween  a  first  guess  field  (in  this  case  the  previous  anal- 
ysis weighted  towards  climatology  at  a  specied  rate)  and  the 
assembly  field,  are  eliminated  in  a  weighted,  least  mean 
squares  sense.  The  first  guess  field  provides  shape  informa- 
tion for  the  horizontal  blending  procedure.  The  assembly 
field,  into  which  all  new  observations  are  assembled,  is 
obtained  in  a  two-step  process.  First,  the  assembly  field 
from  the  previous  analysis  is  trended  towards  climatology  in 
the  same  manner  as  tte  first  guess  field,  and  a  weight  field 
associated  with  it  is  decayed  an  arbitrary  amount.  Then,  a 
weighted  average  of   the  assembly  and  first   guess  fields  is 
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The  vertical  grid  on  which  the  twenty-six  EOTS 
parameters  are  defined.  Columns  labeled  T,  T',  T"  are 
temperature,  first  vertical  temperature  derivative,  and 
second  vertical  temperature  derivative,  respectively. 
Parameters  y-2b  are  associated  with  fixed  levels. 
Parameters  2-8  are  associated  with  floating  levels 
(PLD)^  relative  to  Parameter  1,  Primary  Layer  Depth 


Figure  2.1    Analysis  Levels  for  FNOC  Analyses. 
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formed  (using  a  low,  specified  weight  for  the  first  guess 
field)  to  provide  a  new  assembly  field.  New  observations  are 
interpolated  to  the  nearest  analysis  grid  point  using  shape 
information  from  the  first  guess  field.  At  each  grid  point, 
a  weighted  average  of  the  interpolated  observations  and  the 
assembly  field  is  formed  and  a  final  assembly  field  deter- 
mined. It  is  upon  this  final  assembly  field  and  the  first 
guess  field  that  the  PIB  analysis  is  performed.  The  proce- 
dure is  designed  to  spread  information  from  data-rich  points 
to  data-sparse  points  and  produces  a  smoothed,  horizontally 
and  vertically  blended  field.  During  the  blending,  both 
first  guess  and  assembly  fields  are  weighted  towards  clima- 
tology, and,  in  the  absence  of  inputs,  they  will  revert  to 
climatclcgy  over  a  period  of  approximately  fifteen  days 
[Pollack,  personal  communication]. 

The  analysis  is  performed  on  twenty-six  parameters 
at  fixed  and  floating  levels  from  the  surface  to  400m,  Fig 
2. 1  .  Parameter  1  is  the  Primary  Layer  Depth  (PLD)  which  in 
general  eguates  to  the  depth  of  the  surface  mixed  layer.  The 
remaining  parameters  are  temperature  at  selected  depths  and 
vertical  temperature  gradients.  Parameters  2  to  8  are  anal- 
ysed at  floating  levels  defined  relative  to  the  PLD. 
Parameters  9  to  26  are  associated  with  fixed  levels.  Below 
400m  the  thermal  field  is  constructed  from  climatological 
archives  modified  to  blend  smoothly  with  the  analysed  temp- 
erature profiles  above. 

Selection  of  the  PLD  (point  1)  is  accomplished  by 
applying  an  algorithm  to  all  BT  observations  made  during  a 
specified  number  of  days  (typically  five)  prior  to  the  anal- 
ysis time  to  identify  PLD  candidates  (i.e.,  the  depth  (s)  at 
which  the  criteria  prescribed  by  the  PLD  algorithm  are 
satisfied)  from  each  profile.  The  algorithm  assigns  the  PLD 
to  the  depth (s)  (below  the  surface)  of  maximum  rate  of 
curvature  in   the   vertical  profile.    Then  a   preliminary 
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three-cycle  FIB  is  performed  using  the  selected  PLD 
candidates  as  input,  and  the  PLD  and  assembly  fields  from 
the  previous  analysis,  trended  towards  climatology,  as  the 
first  guess  and  initial  assembly  fields.  Any  recent  BT's  are 
then  examined  and  the  PLD  candidate  from  each  that  is 
closest  to  the  PLD  produced  by  the  preliminary  FIB  analysis 
is  selected  and  entered  into  a  second  three-cycle  FIB 
analysis.  This  produces  the  final  analysed  PLD. 

2.   TOPSrCoupled  FOTS  (TEOTS)  Analysis 

The  TEOTS  system  is  an  attempt  to  incorporate 
oceanic  physics  and  air-sea  interactive  processes  into  the 
analysis  by  coupling  the  ocean  thermal  analysis  to  the 
atmospheric  forcing  via  the  physics  of  the  TOPS  mixed  layer 
model. 

She  TOPS  mixed  layer  model  is  the  Navy's  vertical 
ocean  thermal  structure  forecast  system.  Given  that  the 
ocean  is  primarily  atmospherically  forced,  and  that  mixed 
layer  depths  have  proved  'highly  predictable  with  a  variety 
of  models',  such  as  the  Denman  (1973)  or  Mellor  and  Durbin 
(1975)  models,  TOPS  was  developed  by  NORDA  (Naval  Ocean 
Research  and  Development  Activity) ;  it  is  based  largely  on 
the  use  of  rather  well-defined  atmospheric  variables  to 
drive  a  model  to  describe  the  relatively  data  sparse  oceans. 
Surface  wind,  solar  heat  fluxes,  and  precipitation  fields 
are  among  those  necessary  to  drive  the  TOPS/TEOTS  models. 
The  surface  fluxes  are  supplied  by  NOGAPS  (Naval  Operational 
Global  Atmospheric  Prediction  System)  and  the  physics  of  the 
model  are  based  on  the  Mellor  and  Yamada  (1974)  level-2 
turbulence  model  with  the  added  effects  of  advection  by 
instantaneous  wind  drift  and  climatologically  averaged  geos- 
trophic  currents.  The  TOPS  model  produces  a  twenty-four 
hour  forecast  for  temperature  fields  at  seventeen  depths 
from   the   surface  to  500m   on   a  horizontal   grid   spacing 
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identical  to  that  used  in  the  teots  analysis  to  which  it  is 
coupled-  The  main  purpose  behind  the  TOPS/TEOTS  approach  is 
to  de-couple  the  objective  analysis  component  (EOTS)  SST  and 
MLD  fields  from  climatology,  and  to  provide  realistic  lower 
boundary  conditions  fcr  the  atmospheric  models. 

TEOTS  is  coupled  to  TOPS  in  cyclical  fashion, 
providing  initial  conditions  for  a  24hr  TOPS  forecast  that 
is  subsequently  fed  back  into  EOTS  as  a  first  guess  field 
for  the  following  day's  analysis.  With  the  exception  of  the 
lock-step  procedure  described  above  and  a  different 
prescription  for  certain  parameters,  the  TEOTS  analysis 
system  is  the  same  as  EOTS,  relying  on  the  FIB  technique  to 
combine  the  various  types  of  data.  Clancy  and  Pollack  (1983) 
give  a  mere  rigorous  description  of  both  EOTS  and  TEOTS. 

For  the  study  domain  EOTS  was  used  for  the  regional 
(NOCAL)  analysis  and  TEOTS  for  the  hemispherical  analysis. 

B.   ACOUSTIC  HODE1S 

1-  lk£   ZACT 9 H  Mcdel 

The  FACT9H  Transmission  Loss  Model  is  the  Navy 
Interim  Standard  Transmission  Model  for  ocean  regions  which 
may  be  treated  with  a  single  sound-speed  profile  and  a  flat 
ocean  bottom  [Jacobs,  1982].  It  is  a  ray-acoustics  model 
designed  for  the  computation  of  transmission  loss  as  a  func- 
tion cf  range  and  frequency  for  a  source  and  receiver  at 
fixed  depths.  The  classical  ray  treatment  has  been  augmented 
with  the  addition  of  higher  order  asymptotic  corrections  in 
the  vicinity  of  caustics,  and  the  phase  addition  of  certain 
ray  paths. 

Wave  theory  for  acoustic  propagation  states  that 
there  is  a  frequency  below  which  no  'trapped'  nornal  mode 
may  exist  for  a  given  channel  or  duct  (waveguide)  [Coppens 
et  al.,   1980],   The  FACT9H  model  attempts  to  translate  this 
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condition  to  ray-acoustics  by  applying  an  algorithm  [Clay, 
1968  ]  to  characterize  the  gross  features  of  ducted  propaga- 
tion. The  intensity  in  the  surface  duct  is  found  from  the 
principle  of  conservation  of  energy  modified  by  additional 
losses  (proportional  to  range)  caused  by  duct  leakage  and 
rough-surface  scattering  of  energy  from  the  duct  [Marsh  and 
Schulkin,  1967].  For  both  source  and  receiver  within  the 
surface  duct  the  transmission  loss  (TL)  as  a  function  of 
range  (R)  is  given  by  Egn.  2. 1 

TL(R)=22  +  10log(R.9)  ♦  b.R  (egn  2.1) 

where  e  is  the  angle  for  the  ray  at  the  surface  which  just 
grazes  the  bottom  of  the  duct  and  b  contains  the  duct- 
leakage    and  rough-surface   losses, Egn.    2.2 

b=14.88  x    105    f-5/3    g-i/3    z-3    ♦   fl   f V2    z ~ V2  (egn   2.2) 

In  the  above  expression  f  is  the  frequency,  g  the  magnitude 
of  the  below  layer  gradient,  z  the  duct  depth,  and  W  is  a 
factor  determined  by  surface  wave  height. 

The  9H  is  a  model  designator  and  reflects  the  stage 
of  development  of  the  model  physics  (FACT10  is  presently  in 
operation  at  FNOC) . 

2.   The  PE  Model 

The  PE  Model  replaces  the  reduced  elliptic  wave 
equation  with  a  parafcclic  partial  differential  equation  that 
can  be  integrated  numerically  using  the  Tappert-Hardin 
split-step  Fourier  Algorithm  [Brock,  1978].  The  parabolic 
wave  equation  includes  diffraction  and  all  other  full-wave 
effects  as  well  as  range  depedent  environments.  The  entire 
range  and  depth  dependent  acoustic  field  is  computed  as  the 
solution  is  marched   forward  in  range.   The   model  assumes  a 
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flat  pressure  release  ocean  surface  and  a  vanishing  field  at 
the  depth  of  the  finite  Fourier  transform  (i.e.,  for  bottom 
grazing  angles  greater  than  the  maximum  permissible  33° ,  the 
acoustic  energy  is  smoothly  attenuated  to  zero  and  back 
reflection  of  acoustic  energy  into  the  water  column  is 
avoided) .  For  grazing  angles  less  than  the  maximum,  a 
pseudo  radiation  condition  is  introduced  at  the  water-bottom 
interface  by  smoothly  attenuating  the  field.  Since  the  error 
in  the  parabolic  approximation  increases  as  angles  increase 
from  the  horizontal,  steep  bottom  slopes  can  cause  inaccura- 
cies. The  model  is  primarily  considered  useful  for 
predicting  low- frequency  acoustic  propagation  of  energy 
along  waterborne  or  shallow-angle,  bottom-bounce  paths  in 
deep  water. 
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III.    DATA    BASE 

A.  DATA  ACQUISITION 

The  XBT  and  CTD  data  from  the  surface  to  500m  acquired 
from  the  OPTOMA5  cruises  formed  the  primary  source  of  field 
data.  CPTCMA5  consisted  of  three  separate  cruises  undertaken 
between  15  June  and  19  July,  1983  by  R/  V  ACANIA.  The  three 
cruises,  designated  AI,  All  and  AIII  each  sampled  approxi- 
mately the  same  area  of  the  CCS,  Fig.  3. 1 

A  secondary  source  of  field  data  was  provided  by  deep 
CTD  casts  from  the  cruise  undertaken  by  the  R/V  DE  STEIGUER 
between  10  and  21  June,  1983.  A  section  of  this  cruise 
sampled  the  0PT0MA5  area  and  provided  CTD  data  to  3000m. 
Satellite  imagery  was  also  used,  when  available,  to  provide 
a  qualitative  measure  of  the  sea  surface  temperature 
distribution. 

Temperature  profiles,  to  be  compared  later  with  the 
field  data,  were  retrieved  from  the  archives  of  FNOC  for  the 
same  area  and  period.  Both  EOTS  and  TEOTS  depth/temperature 
profiles  for  the  beginning  and  end  of  each  OPTOMA5  cruise 
were  obtained  for  selected  positions,  Fig.  3.2  . 

Details  of  the  data  base  are  given  in  Table  I  . 

B.  DESCRIPTION  OF  DATA  ABO  DATA  PROCESSING 

1.   ACANIA  Data 

The  positions  of  the  XBT/CTD  casts  for  the  Acania 
cruises  are  documented.  Figs.  3.3  to  3.5.  All  XBT/CTD  data 
were  edited  thus  eliminating  any  obviously  erroneous 
profiles.  A  secondary  visual  inspection  of  each  profile  was 
undertaken  with  the  result  that  a  final  data  base  for  each 
cruise  was  obtained,  Table  I 
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TABLE  I 
Details  of  Data  Base 


SOUECE 

PEEIOD 

NO.    OF    ACCEPTABLE 
PROFILES 

RANGE   OF 
PROFILES 

OPTOMA    AI 

6/15-6/22 

74 

0-500m 

OPTOMA    All 

6/29-7/4 

89 

0-50  0m 

OPTOMA    AIII 

7/13-7/19 

86 

0-50  0m 

DE    STEIGUER 

6/10-6/21 

8 

0-30  00  m 

EOTS/TEOTS 

6/16,6/23, 

6/30,7/7. 

7/10,7/21 

42 

0-40 0m 

The  XBT's  were  then  interpolated  in  the  vertical  to 
standard  5m  depths  from  the  surface  to  500m.  The  loss  of 
temperature  microstructure  resulting  from  this  interpolation 
was  assumed  unimportant  in  view  of  the  goals  of  this  study. 
These  edited,  'smoothed'  depth/temperature  profiles  formed 
the  field  data  upon  which  the  later  analysis  was  to  be 
based . 

Information  from  each  cruise  was  assumed  'quasi- 
synoptic'  (i.e.  the  region  was  believed  to  have  varied  some- 
what over  the  period  of  sampling) .  This  assumption  is 
weakest  in  the  near  surface  layers  where  the  time  scale  of 
variability  can  be  of  the  order  of  a  day  or  less.  The 
validity  of  this  assumption  in  the  context  of  the  present 
study  is  addressed  in  Chapter  V. 

2.   EZ  STEIGDER  Data 

The  DE  STEIGDER  sampling  provided  CTD  casts  from  the 
surface  to   a  maximum  depth   of  3000m.   These  profiles  were 
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also  interpolated  vertically  to  5m  depth  intervals  and, 
based  on  the  assumption  that  the  water  column  was  near 
stationary  below  the  therraocline  (approximately  500m) ,  they 
were  used  to  extend  the  ACANIA  and  FNOC  profiles  to  a  depth 
of  3000m.  These  profiles  were  then  extrapolated  linearly  to 
the  ocean  bottom  (assumed  flat  at  4000m) . 

An  average  salinity/depth  profile  was  calculated 
from  the  'in-area'  EI  STEIGUER  CTD  casts  and  this  average 
profile  was  assumed  as  the  standard  salinity  profile  for  all 
sound  speed  calculations. 

3 .   Satellite  Imagery 

Advanced  Very  High-Resolution  Radiometer  (AVHRR) 
infrared  images  from  the  NOAA-7  satellite  were  used  to  qual- 
itatively show  the  icre  obvious  features  of  the  SSI  field. 
These  images  were  obtained  on  an  'as  available1  basis  during 
each  of  the  ACANIA  cruises.  Details  of  the  major  features 
shown  in  these  images  were  sketched,  Figs.  3.6  to  3.8  The 
satellite  images  were  also  used  to  gualitively  assess  the 
assumption  of  guasi-synopticity. 

U.   FNOC  Archived  Data 

Temperature  profiles  for  the  standard  positions, 
Fig.  3.2  ,  were  extracted  from  the  archives  of  FNOC.  These 
profiles  were  used  fcr  both  EOTS  and  TEOTS  analyses  on  the 
dates  closest  to  the  commencement  and  termination  of  each 
ACANIA  cruise. 
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IV.  DATA  ANALYSIS 

A-   EATA  PBOCESSIHG 

To  facilitate  ccnparison  with  FNOC's  data  fields,  the 
study  domain  was  divided  into  a  7x7,  square  grid  with  a 
grid-spacing  of  20km  and  the  irregularly  spaced  XBT/CTD 
profiles  interpolated  onto  this  grid. 

The  horizontal  correlation  of  the  temperature  field  at 
various  depth  intervals  in  the  top  150m  of  the  water  column 
were  evaluated.  The  correlation  at  various  levels  displayed 
similar  chracteristics  showing  zero  crossings  occurring  at 
approximately  80km  with  a  0.5  correlation  at  approximately 
30km  (approximately  the  e-folding  distance  between  the 
origin  and  the  zero  crossing  distance).  It  was,  therefore, 
reasonable  to  assume  that  the  temperature  variations  at 
points  within  30  km  cf  a  given  position  within  the  domain 
were  related  (i.e.  all  points  within  30km  could  be  expected 
to  influence  a  given  point) .  As  an  example,  the  SST  hori- 
zontal correlations  for  each  cruise  are  presented,  Figs.  4.1 
to  4.3  .  The  correlation  information  thus  derived  was  used 
when  interpolating  the  field  data  onto  the  grid.  Only 
observations  within  30km  of  each  grid  point  were  used  and 
each  observation  value  was  weighted  according  to  the  inverse 
of  its  distance  from  the  grid  point.  The  interpolated  data 
thus  formed  were  then  automatically  contoured  at  FNOC  temp- 
erature levels,  Fig.  2.1  The  mixed  layer  depth  (MLD)  calcu- 
lated frcm  the  field  data  was  equated  to  the  PLD  in  the  FNOC 
analyses  (point  1  Fig.  2.1)  and  was  defined  as  the  shal- 
lowest depth  below  which  a  negative  temperature  gradient  was 
found  with  a  magnitude  greater  than  3°C  /100m  over  two 
consecutive  5m   depth  intervals.   The   algorithm   used   to 
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calculate  the  MLD  was  designed  to  yeild  the  true  MLD  even  in 
the  presence  of  strong,  near-surface  negative  temperature 
gradients  caused  by  diurnal  heating  the  "afternoon  effect" 
which  seldom  exceeded  7m  in  depth.  MLD's  calculated  by  this 
method  agreed  well  with  manually  determined  MLD's  for  the 
sampled  profiles.  PNOC  EOTS  and  TEOTS  profiles  were  also 
interpolated  and  contoured  in  a  similar  manner  and  on  the 
same  grid.  The  MIE  computed  for  the  model  profiles  using 
the  algorithm  described  above  was  in  agreement  with  the  PLD. 

By  contouring  the  fields  in  the  manner  described,  it  was 
possible  to  extract  a  depth/temperature  profile  from  the 
surface  to  4030m  for  any  position  within  the  study  domain. 

Throughout  the  study  domain,  temperature  variations  at 
400m  were  small  enough  (<0.5°C)  for  all  profiles  (including 
those  provided  by  FNOC)  to  be  smoothly  blended  between 
shallow  and  deep  profiles  above  the  depth  of  the  deep  sound 
channel  axis  (ca.  500m)  . 

The  feasibility  of  using  discrete  depth/temperature 
levels,  with  depth  separations  dictated  by  the  EOTS  analysis 
levels,  to  represent  the  vertical  temperature  structure  of 
the  area  was  tested  by  examining  the  vertical  pattern  corre- 
lation between  temperatures  in  the  near  surface  layers  (the 
zone  of  maximum  variation  of  temperature  with  depth) ,  Figs. 
4.4  to  4.6  .  These  pattern  correlations  are  similar  for  each 
cruise  and  have  zerc  crossings  at  depths  of  the  order  of 
300m  when  the  SST  is  correlated  with  deeper  levels.  Such  a 
pattern  correlation  gives  an  good  indication  of  the  prob- 
ability that  variations  in  temperature  at  one  depth  are 
related  to  variations  in  temperature  at  another  depth.  At 
the  e-folding  depth  (approximately  100m)  ,  the  correlation  is 
approximately  0.7.  It  was  therefore  considered  reasonable  to 
represent  the  vertical  temperature  structure  at  the  discrete 
levels  dictated  by  the  EOTS  analysis.  The  EOTS  depth 
interval   between  temperature   levels   in   the  near   surface 
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Figure  4. 1    SSI   Horizontal  Correlation-Cruise  AI 
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Figure   a. 2        SST   Horizontal  Correlation-Cruise  All. 
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Figure   4.3        SST   Horizontal  Correlation-Cruise  AIII- 

U3 


layer  is  of  the  order  of  50m  which  yields  a  pattern  correla- 
tion coefficient  ranging  from  0.6  for  cruise  All  to  0.8  for 
cruises  AI  and  AIII .  Below  the  mixed  layer  (approximately 
50m  )  the  correlation  improves  to  0.8  or  better  between 
adjacent  contoured  levels. 

The  smoothing  effect  of  deriving  profiles  as  described 
above  (i.e.  the  higher  frequency  fluctuations  of  temperature 
with  depth  are  removed)  was  not  considered  important, 
because  the  profiles  thus  derived  were  later  to  be  used  as 
inputs  tc  low  frequency  (<1000Hz)  acoustic  models  and  the 
essential  shape  cf  the  profile  was  maintained. 
Additionally,  ray  theory  is  only  valid  if  sound  velocity 
does  not  change  much  in  a  wavelength,  therefore,  the  inclu- 
sion of  the  higher  wavenumber  fluctuations  of  temperature 
(and,  therefore,  sound  speed)  with  depth  could  invalidate 
the  use  of  this  theory  [Urick,  1983]. 

Four  features  of  the  vertical  temperature  profile  were 
selected  as  being  of  primary  importance  in  describing  the 
thermal  structure  for  low  frequency  passive  acoustic  appli- 
cations. These  were:  sea  surface  temperature  (SST)  ,  mixed 
layer  depth  (MLD)  (previously  defined) ,  the  below  layer 
gradient  (BLG)  (defined  as  the  temperature  gradient  in  the 
25m  immediately  below  the  MLD)  and  the  thermocline  gradient 
(THG)  (Defined  as  the  temperature  gradient  between  the  MLD 
and  the  deep  sound  channel  axis  depth) .  The  contour  fields 
derived  for  each  of  these  parameters,  Figs.  4.7  to  4.10  , 
are  described  in  detail  in  Section  B  of  this  chapter. 

B.   DESCRIPTION  OF  ARIA  THERMAL  STRUCTURE 

The  following  is  a  brief  description  of  the  thermal 
structures  measured  in  the  study  domain  during  the  three 
CPT0MA5  cruises.  The  desciption  in  limited  to  the  four 
parameters  defined  in  section  A  above. 
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Figure  4.4   Pattern  Correlation-Cruise  AI 
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Figure   4.5        Pattern   Correlation-Cruise   All. 
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Figure  4.6        Pattern   Correlation-Cruise  AIII 
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a. 


b. 


c. 


Contour  units:      C  x   10 


Figure   4.7        SST  Contours 
a.    Cruise   Ai      b.    Cruise  All      c.    Cruise  AIII. 
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Contour  units:  metres 


Figure  4.8    HLD  Contours 
a.  Cruise  hi     b-  Cruise  All   c.  Cruise  AIII. 
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Figure  4.9    BLG  Contours 
a.  Cruise  At   b.  Cruise  All   c.  Cruise  AIII. 
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Figure  4. 10    THG  Contours 
a.  Cruise  Ax   b.  Cruise  All   c.  Cruise  AIII. 
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The  significant  weather  for  each  cruise  (i.e.  wind  speed 
and  clcud  cover)  were  obtained  from  R/V  ACANIA's  hourly  log. 
The  wind  was  from  a  generally  southeastward  throughout  all 
three  cruises,  and,  for  the  most  part,  wind  speeds  were 
between  15  and  20  knots.  The  major  excursions  from  this  norm 
were  (a)  at  the  beginning  of  cruise  AI  when  wind  speeds  were 
greater  than  30  knots  on  the  first  and  second  days  (16  and 
17  Jun),  (b)  during  the  latter  part  of  cruise  AI  and  early 
in  cruise  All  when  winds  of  10  to  15  knots  were  recorded  and 
(c)  during  the  first  two  days  of  cruise  AIII  when  wind  speeds 
again  increased  to  30  knots  or  more.  Throughout  the  three 
cruises,  partly  cloudy  to  clear  skies  prevailed,  with 
predominantly  less  than  4  oktas  of  cloud  cover.  Pressure 
remained  nearly  steady  at  approximately  1030  mb. 

1-   Cruise  A I  Q5  Jun-22  Jun,  1983) 

SST's  ranged  from  14.6°C  in  the  northwest  of  the 
domain  tc  12.2°C  associated  with  a  cool  feature  protruding 
into  the  southeast,  Fig.  4.7  (a) .  The  effect  of  the  cool 
intrusion  was  to  intensify  the  horizontal  temperature 
gradient  in  its  vicinity  and  to,  locally,  distort  the  other- 
wise nearly  meridional  SST  structure.  The  maximun  tempera- 
ture gradient  intensification  occurred  on  the  northern  flank 
of  the  cool  feature,  where  a  gradient  of  7.0oC/100n  was 
observed.  The  position  of  the  cooler  water  agreed  well  with 
the  position  obtained  from  satellite  imagery.  Fig  3.6. 

The  depth  of  the  mixed  layer.  Fig.  4.8  (a)  ,  varied 
from  20  to  55m  in  a  complex  manner  throughout  the  area  with 
the  shallowest  depths  in  a  50km-wide  swath  southeast  to 
northwest  across  the  centre  of  the  area.  No  discernable 
correlation  was  observed  between  the  SST  and  MLD  patterns. 

Below  layer  gradients  ranged  from  3.5°C/100m  in  the 
sector  affected  by  the  cold  feature  in  the  southeast  to 
5.0°C/100m   in   an   area   apparently  associated   with   the 
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shallower  MLD's,  Fig  4.9  (a).  Thermocline  gradients  varied 
in  a  near  meridional  manner  with  the  weaker  gradients  in  the 
west.  The  magnitude  of  the  gradient  varied  from  1.5°C/100m 
to  1.8°C/100m,  Fig  4.10  (a). 

2.   Cruise  All  (29  Jiin-4  Jul,  1983) 

A  general  warming  trend  was  apparent  between  the  end 
of  cruise  AI  and  the  beginning  of  cruise  All  with  SSI's 
generally  2.0°C  warmer  throughout  the  domain.  There  was  also 
a  distinct  change  in  the  SST  pattern.  Fig  4.7  (b) .  The  cool 
feature  initially  observed  in  the  southeast  sector  had 
relaxed  southwards  and  a  second  cool  feature  was  observed  to 
extend  northward  into  the  southwest  of  the  domain.  The 
warmest  water  remained  in  the  far  northwest  sector,  with  a 
near  meridional  structure  throughout  the  remainder  of  the 
area.  Temperatures  ranged  from  14.4°C  to  16.2°C.  The  warming 
trend  was  to  be  expected  for  the  period  involved. 
Climatology  [Podeszoa,  1976]  shows  the  SST  increasing  from 
14.1°C  in  June  to  16.5°c  in  July.  The  mainly  light  to  fresh 
winds  and  clear  skies  prevailing  through  the  latter  part  of 
AI  and  early  days  of  All  could  be  expected  to  produce  the 
warming  trend  noted.  Satellite  imagery  once  more  agreed 
well,  gualititavely,  with  the  position  of  the  major  surface 
features  shown  by  the  analysis.  Fig  3.7. 

Associated  with  the  warming  trend  and  generally 
lighter  winds,  the  HID  shoaled  over  the  period  between  AI 
and  All.  Layer  depths,  some  20m  shallower  than  than  those 
observed  during  AI,  were  measured  with  depths  varying  from 
near  surface  in  the  northwest  to  40m  in  the  southeast  but 
with  a  complex  structure  throughout  the  area. 

The  strengthening  of  the  BLG,  Fig  4.9  (b)  compared 
to  those  observed  during  cruise  AI  indicated  that  the 
warming  trend  had  been  limited  to  the  upper  layers  of  the 
ocean  and  was,   therefore,   probably  atmospherically  driven. 
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Magnitudes  ranged  frcm  10. 0oc/100m  in  the  west  in  associa- 
tion with  the  lower  SST's  to  3.5°C/100m  in  the  centre  of  the 
domain,  loosely  associated  with  the  deeper  MLD's.  An 
increase  in  THG's  to  tetween  1.7oc/100m  and  2.1°C/100m  with 
a  generally  similar  pattern  to  that  of  cruise  AI,  emphasized 
the  supposition  that  the  warming  had  been  limited  to  the 
upper  layers  of  the  ocean. 

3-   Cruise  AIII  (J3  Jul^lS  Jul,  1283) 

The  SST  analysis,  Fig.  4.7  (c)  showed  the  presence 
of  cool  water  along  the  southern  edge  of  the  domain  with  the 
cool  feature  in  the  southwest  during  cruise  All,  having 
relaxed  southward  and  a  further  intrusion  of  cocl  water 
extending  northward  into  the  area  from  the  southeast. 
Satellite  imagery  agreed  well  with  the  position  of  these 
features,  Fig.  3.8  .  However,  coverage  was  too  intermittent 
to  allow  the  feature  in  the  southeast  to  be  positively 
related  to  the  cool  feature  measured  in  a  similar  position 
during  cruise  AI .  Temperatures  ranged  from  12.4  to  16.0°C 
with  the  warmest  water  remaining  as  an  almost  stationary 
feature  in  the  far  northwest  of  the  domain.  A  relatively 
intense  horizontal  temperature  gradient  of  10.0°C/100m  was 
measured  in  a  30km-wide  band  around  the  cool  feature  in  the 
southeast. 

MLD's,  Fig.  4.8  (c)  were,  on  average,  deeper  than 
those  measured  during  cruise  All,  indicative  of  the  higher 
winds  prevailing  during  the  early  part  of  the  cruise 
producing  significant  mixing  in  the  surface  layers.  However, 
shoaling  of  the  MLD,  in  both  the  north  and  south  of  the 
area,  tc  10m  or  less  would  tend  to  indicate  that  the  wind 
mixing  was  restricted  to  the  area  being  sampled  by  ACANIA 
during  that  period.  The  general  horizontal  structure  of  the 
MLD  remained  complex. 
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The  magnitude  of  neither  the  BLG  nor  the  THG  changed 
much  from  those  measured  during  cruise  All.  However,  the 
distribution  of  the  iraximum  and  minimum  BLG  was  signifi- 
cantly different  from  that  of  All,  although  there  was  a  weak 
increase  in  BLG  in  association  with  the  cool  feature  in  the 
southeast  of  the  domain. 
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V.  COMPARISON  OF  OBSERVED  AND  MODELED  TEMPERATURE  STRUCTURE 

For  each  of  the  parameters  discussed  in  the  Chapter  IV, 
namely,  SST,  MLD,  BLG  and  THG,  a  difference  field  was 
produced  by  subtracting  the  magnitude  of  the  parameter 
derived  from  the  model  analysis  from  that  derived  from  the 
0PT0MA5  data.  The  procedure  was  carried  out  for  each  model 
analysis  at  each  of  the  49  grid  points  in  the  domain.  The 
difference  fields  thus  calculated  were  contoured,  Figs.  5.1 
to  5.24. 

Statistics  were  calculated  from  the  contoured  difference 
fields  by  examining  the  percentage  of  the  total  domain 
enclosed  by  discrete  contour  intervals  in  each  difference 
field,  and  for  each  parameter.  The  • bin-size'  for  each  field 
was  established  by  examining  the  individual  plots  and 
subjectively  selecting  an  appropriate  difference  interval 
such  that  each  bin  would  represent  between  10  and  20*  of  the 
dynamic  range  of  the  variable,  whilst  keeping  the  area 
covered  by  each  bin  large  enough  to  make  measurement  prac- 
tical. The  bin-sizes  selected  were: 

SST  -  0.4°C 

MLD  -  10.0m 

BLG  -  1.0OC/100m 

THG  -  0.  1oc/100m 
The  areas  within  each  bin   were  measured  using  a  Planimeter. 
Figs.   5.25  to  5.48  show  the  histograms  of  the  distributions 
obtained  by  the  above  method  and   Tables  II  to  IV   list  the 
relevant  statistics  for  each  histogram. 

A  negative  value  in  the  difference  calculations  implies 
that  the  model  analysis  was  overestimating  that  particular 
parameter  (i.e.,  Difference  =  Sea  truth  -  Model  value). 


56 


Because  of  a  corruption  in  the  data  supplied  by  FNOC  it 
was  not  possible  to  derive  contour  fields  of  the  model  anal- 
yses for  the  termination  of  cruise  AIII. 

The  larger  grid  spacing  used  in  the  TEOTS,  hemispherical 
analyses  (ca.  320km)  meant  that  TEOTS  comparisons  were, 
necessarily,  general  in  nature  and  should  be  viewed  in  terms 
of  average  values  ,  since  it  is  unlikely  that  such  a  large 
grid  spacing  will  allow  the  analysis  to  retrieve  the  mesos- 
cale  variability  present  in  the  area  (  horizontal  vari- 
ability on  a  scale  of  the  order  of  tens  of  kilometers) . 

The  following  sections  discuss  the  statistics  obtained 
and  qualitatively  estimate  the  effectiveness  of  the  model 
analysis  in  representing  the  sea-truth.  A  summary  of  the 
model  performances  is  given  in  Section  A. 6  of  this  chapter. 
Throughout  the  Chapter,  sea  truth  refers  to  0PT0MA5  derived 
data. 

A.   COMPARISONS 

1 .   Cruise  A I 

The  analysis  dates  for  cruise  AI  were  16  June,  1983 
(initial)  and  23  June,  1983  (final) .  FNOC  data  files  showed 
inputs  into  the  the  analyses  for  the  study  domain  of  three 
BT's  for  the  initial,  and  twenty-three  ET's  (from  R/V 
ACANIA)  for  the  final  analyses.  The  statistics  for  the 
difference  fields  derived  for  this  cruise  are  given  in  Table 
II  and  the  histograms  from  which  these  statistics  were 
derived  are  shown  in  Figs.  5.25  to  5.32 

a.   EOTS  Analysis 

The  initial  SST  analysis,  Fig.  5.1  (b)  showed  a 
nearly  meridional  structure  with  the  warmest  water  in  the 
west  of  the  domain  and  a  weakly  strengthening  horizontal 
temperature  gradient  in  the  east.  The  difference  field,  Fig. 
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5.1  (c)  showed  the  analysis  to  be  overestimating  the  meas- 
ured SST  for  the  majority  of  the  area.  The  mean  difference 
was  -0.8°C  and  the  maximum  difference  was  -2.0<>C  associated 
with  the  intrusion  cf  cool  water  into  the  southeastern 
sector.  Fig.  5.1  (a).  The  final  analysis,  Fig.  5.1  (d) 
showed  a  general  movement  of  the  analysis  field  towards  the 
sea  truth  in  both  shape  and  magnitude.  The  shape  field  was 
only  marginaly  improved  and  continued  to  neglect  the  pres- 
ence of  the  cool  feature  in  the  southeast.  The  mean  differ- 
ence reduced  to  0.1 °C  and  the  spread  of  the  difference  field 
about  this  mean  was  also  reduced  compared  to  that  of  the 
initial  analysis.  The  maximum  difference  remained  in  the 
vicinity  of  the  cool  feature  but  with  a  reduced  magnitude  of 
-0.80C. 

The  MLD  was  overestimated  in  both  initial  and 
final  analyses,  Fig  5.2  (b)  and  (d)  .  The  slight  shoaling  of 
the  layer  between  the  two  analyses  (2m)  had  no  significant 
effect  on  the  difference  fields.  Fig.  5.2  (c)  and  (e) .  The 
modeled  MLD  was,  on  average,  18  to  19.5m  greater  than  the 
sea  truth  and  the  maximum  difference  was  between  30  and  40m. 
The  near  uniform  structure  of  the  modeled  MLD  in  both  anal- 
yses was  a  gross  misrepresentation  of  the  complex  structure 
evident  in  the  sea  truth  contour  field,  Fig.  5.2  (a). 

The  average  magnitude  of  both  the  3LG,  Fig.  5.3 
and  1HG,  Fig  5.4  model  analysis  fields  also  moved  closer  to 
the  sea  truth  between  the  initial  and  final  model  analyses. 
The  BLG  analyses,  however,  showed  a  marked  difference 
between  modeled  and  sea  truth  shape  fields  for  both 
analyses. 

t.   TEOTS  Analysis 

The  initial  TEOTS  SST  analysis.  Fig.  5.5  (b) 
overestimated  the  sea  truth  in  a  similar  manner  to  that 
observed   in    the   eguivalent   EOTS   analysis    and   by 
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approximately  the  same  average  amount  i.e. ,  a  mean 
difference  of  -0.9°C.  The  mean  difference  was  almost  halved 
in  the  final  analysis,  Fig.  5.5  (c)  .  However,  the  error  in 
the  modeled  shape  field  remained  constant  as,  therefore,  did 
the  difference  distribution  about  the  mean.  The  difference 
fields  are  shown  in  Fig.  5.5  (d)  and  (e) . 

The  modeled  MLD  was  unresponsive  to  input  data 
and  showed  a  lack  of  features.  The  TEOTS  analyses  overesti- 
mated the  MLD  by  an  average  of  12m.  TEOTS  MLD's  were  closer 
to  the  mean  sea  truth  MLD  than  were  the  EOTS  analyses  but, 
because  of  the  lack  cf  shape  agreement  between  model  and  sea 
truth  the  distribution  of  the  difference  field  about  the 
mean  was  similar  to  that  observed  in  the  EOTS  analyses. 

The  BLG,  Fig.  5.7  and  THG,  Fig.  5.8  appear  to 
respond  to  the  input  data  by  moving  closer  to  the  sea  truth 
in  the  mean.  However,  the  lack  of  shape  agreement  between 
model  and  sea  truth  contour  fields  remained  unchanged.  The 
TEOTS  analysis  was  initially  further  from  the  sea  truth  and 
responded  more  slowly  to  input  data  than  the  EOTS  analysis. 

2-   Cruise  All 

The  analysis  dates  for  cruise  All  were  30  June,  1983 
(initial)  and  7  July,  1983  (final)  .  The  inputs  to  the  anal- 
yses were  3  BT's  and  nil  for  the  initial  and  final  analyses, 
respectively.  Statistics  for  the  model/sea  truth  comparisons 
are  given  in  Table  III  and  the  histograms  from  which  the 
statistics  were  derived  are  shown  in  Figs.   5.33  to  5.40 

a.   EOTS  Analysis 

The  initial  SST  analysis.  Fig.  5.9  (t)  showed 
similar  shape  characteristics  to  those  observed  in  the  23 
June  EOTS  analysis  (cruise  AI)  and  a  general  warming  trend 
of  0.5°C  was  indicated.  Over  the  same  period  the  sea  truth 
analysis,  Fig.  5.9  (a)  indicated  an  increase  in  SST  of  2. 0°C 
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and  a  marked  change  in  the  shape  of  the  contour  field.  The 
initial  difference  field,  Fig.  5.9  (c)  showed  the  analysis 
to  be  an  underestimate  of  the  sea  truth  by  an  average  of 
1.0°C  with  a  maximum  difference  of  1.6°C  and  a  general  lack 
of  agreeient  between  model  and  sea  truth  shape  fields  was 
apparent. 

The  final  model  analysis.  Fig.  5.9  (d) ,  showed  a 
trend  towards  a  meridional  SST  structure  (away  from  the 
zonal  structure  apparent  in  the  sea  truth  field)  with  no 
change  in  the  average  SST.  The  mean  difference  remained 
similar  to  that  calculated  for  the  initial  analysis. 
However,  because  of  the  increased  disagreement  between  the 
model  and  sea  truth  shape  fields,  the  spread  of  the  differ- 
ence distribution  about  the  mean  increased. 

The  FNOC  climatological  depth/temperature 
profile  for  the  June/July  period  in  the  region  of  the  study 
domain  predicts  a  SST  of  14.2°C  which  was  in  close  agreement 
with  the  mean  SST  of  14.4°C  observed  in  the  model  analysis. 

The  initial  MLD  analysis,  Fig.  5.10  (b)  showed 
no  significant  change,  in  shape  or  magnitude  from  that 
observed  in  the  23  June  analysis.  Fig.  5.2  (d) .  However, 
changes  in  the  sea  truth  field,  Fig.  5. 10  (a)  ,  between 
cruises  AI  and  All  meant  that  the  difference  field,  Fig. 
5.10  (c) ,  was  markedly  different  from  that  of  the  23  June 
analysis.  Although  the  final  MLD  analysis,  Fig  5.10  (d)  did 
indicate  a  slight  shoaling  (5m  to  10m) ,  the  MLD  was  overes- 
timated by  the  model  analyses  throughout. 

The  shoaling  of  the  MLD  between  the  initial  and 
final  analyses  noted  above  is  compatible  with  the  model 
analysis  moving  towards  the  FNOC  July  climatology  for  the 
region  which  shows  a  change  in  MLD  from  55m  in  June  to  37m 
in  July. 

The  model  BLG's  were  overestimated  and  the  shape 
of  the  contour  field  in  poor   agreement  with  that  of  the  sea 
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truth.  A  smoothing  of  the  modeled  BLG  field  to  a  near 
constant  value  in  the  final  analysis  did,  however,  reduce 
the  mean  difference  from  the  1.2°C/100ra  observed  in  the 
initial  analysis  to  0.3°C/100m.  The  results  of  the  compari- 
sons are  shown  in  Fig.   5.11- 

Model  and  sea  truth  THG's  were  in  good  agreement 
for  both  analyses,  Fig.  5.12. 

b.   TEOTS  Analysis 

The  warming  trend  noted  in  the  EOTS  analyses  was 
also  observed  in  the  TEOTS  analyses,  Fig.  5.  13  The  mean 
difference  decreased  from  0.4°C  in  the  initial  TEOTS  anal- 
ysis tc  near  zero  in  the  final  analysis.  The  horizontal 
structure  of  the  TEC1S  analyses  remained  meridional  with  a 
weakly  strengthening  horizontal  temperature  gradient  in  the 
east  of  the  domain  compared  to  the  complex  near  zonal  struc- 
ture observed  in  the  sea  truth  SST. 

Modeled  MID's,  Fig.  5.14  remained  unchanged  by 
any  significant  amount  in  both  magnitude  and  structure  from 
those  of  the  TEOTS  analyses  of  cruise  AI.  The  modeled  MLD's 
were  overestimated  by  between  5m  and  35m,  when  compared  to 
the  sea  truth. 

The  BLG,  Fig.  5.15  was  overestimated  by  the 
model  in  both  the  initial  and  final  analyses  with  the 
magnitude  of  the  mean  difference  increasing  by  the  final 
analysis. 

The  THG,  Fig.  5.  16  was  in  reasonably  good  agree- 
ment with  the  sea  truth  throughout. 

3 .   Cr uise  A  III 

As  previously  stated,  model  analyses  were  not  avail- 
able for  the  termination  of  this  cruise.  There  were  no  BT 
inputs  to  either  model  for  the  10  July,  1983  analysis. 
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The  only  significant  change  in  either  model  was  a 
general  increase  in  the  SST  (0. 5°C)  indicated  in  the  EOTS 
analysis. 

The  significant  changes  in  the  difference  fields. 
Figs.  5.17  to  5.24  were  due  to  changes  in  the  sea  truth 
fields  between  cruises  All  and  AIII. 

The  statistics  of  the  difference  fields  and  the 
relevant  histograms  are  shown  in  Table  IV  and  Figs.  5. 41  to 
5.48  As  with  cruise  All,  the  TEOTS  SST  and  MLD  analyses 
were,  on  average,  in  closer  agreement  with  the  sea  truth 
than  the  EOTS  analyses. 

4.   Effectiveness  of   the  EOTS  Analysis  in  Representing 
Injout  Data 

To  estimate  the  accuracy  of  the  EOTS  analysis  in 
responding  to  input  data  the  author's  analysis  program  was 
revised  to  operate  on  a  4x4  grid  with  a  40km  grid  spacing, 
similar  to  that  used  in  the  EOTS  analysis  program.  The 
twenty-three  BT's  input  to  the  EOTS  23  June  analysis  were 
used  as  input.  The  resulting  SST  contour  field  and  that 
produced  by  the  equivalent  EOTS  analysis  are  shown  in  Fig. 
5.49. 

From  Fig.  5.49  it  can  be  seen  that  the  EOTS  analysis 
did  net  respond  accurately  to  the  input  data.  The  accuracy 
of  the  author's  analysis  program  had  previously  been  veri- 
fied by  comparisons  with  manually  contoured  SST  fields  with 
favourable  results.  The  reason  for  the  lack  of  resolution  in 
the  ECTS  analysis  is  not  fully  understood  but  is  probably 
due  to  excess  weighting  on  the  first  guess  field  used  in  the 
EOTS  FIB  analysis.  This  first  guess  field  provides  the  shape 
information  to  which  any  input  data  is  fitted  in  a  least 
squares  sense. 
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5-   Validity  of  Quasi-Synoptic  Assumption 

lack  of  synopticity  in  the  0PT0MA5  data  proved  to  be 
unimportant  since  the  only  inputs,  of  any  significance,  to 
the  FNOC  models  were  the  twenty-seven  BT  submissions  prior 
to  the  23  June  analyses.  The  fact  that  the  BT»s  were  input 
to  the  models  as  one  data  set  made  them  synoptic  to  the 
computer  model  and,  therefore,  a  synoptic  assumption  is 
valid  for  the  sea  truth  analysis  with  which  the  model  anal- 
yses were  compared. 

6 .   Summary 

a.  SST 

The  EOTS  SST  analysis  responded  to  the  magnitude 
changes  dictated  by  input  data  (when  sufficient  data  was 
available)  but  horizontal  structure  was  not  resolved  on  a 
scale  commensurate  with  the  resolution  capability  of  the 
UOkm  grid  spacing  used.  In  the  absence  of  input  data  the 
EOTS  SST  moved  towards  climatology  and  became  even  less 
structured. 

There  was  a  measureable  response  to  input  data 
in  the  TEOTS  SST  analysis.  The  magnitude  of  the  response  was 
less  than  that  observed  in  the  EOTS  analysis  and  the  lack  of 
agreement  in  the  shape  field  between  model  and  sea  truth  was 
even  more  pronounced,  as  might  have  been  expected  frcm  the 
larger  grid  spacing  used  in  the  TEOTS  analysis.  TECTS 
perfomed  better,  in  a  general  sense,  than  EOTS  when  data 
were  not  available  and  the  move  towards  climatology  noted  in 
the  ECTS  analysis  was  not  observed. 

b.  MLD 

Neither  model  produced  a  realistic  analysis  of 
the  sea  truth  MLD  whether  input  data  were  available  or  not. 
Both  analyses   tended  to  overestimate   the  magnitude   of  the 
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MLD,   with  the   TEOTS  analysis  giving  a  consistently  better 
estimate  of  the  sea  truth  MLD  than  the  EOTS  analysis. 

c.   BLG  and  THG 

The  lack  of  agreement  between  modeled  and  sea 
truth  MLD  made  comparisons  of  BLG  and  THG  difficult  (Both 
parameters  were  defined  in  terms  of  the  layer  depth).  Any 
agreement  noted  between  these  two  parameters  was  therefore 
regarded  as  purely  coincidental. 

B.   INTEB-ANHUAL  VARIABILITY 

The  mesoscale  variability  of  the  study  domain  has  been 
dealt  with  in  some  detail  in  the  previous  sections. 
Underlying  this  mesoscale  variability  is  a  variability  with 
a  time  scale  of  the  order  of  years  (the  effect  of  El  Nino) . 
The  region  was  known  to  be  influenced  by  such  an  event 
during  the  period  of  the  present  study.  To  illustrate  the 
effect  of  such  an  anomaly,  mean  depth/temperature  profiles 
were  compared  for  July,  1983  and  July/August,  1982  in  the 
study  domain.  The  mean  profiles  were  also  compared  to  FNOC 
July  climatology  for  the  region.  The  mean  profile  for  1982 
was  derived  from  XB1  data  collected  during  the  0PT0MA2 
cruise  of  July/August,  1982  (The  area  of  the  CCS  sampled  by 
this  cruise  was  very  close  to  that  sampled  in  the  0PT0MA5 
series  of  cruises)  and  that  for  1983  was  derived  from  XBT 
data  collected  during  cruise  AIII  of  the  0PT0MA5  series 
(July,  1983).  Figs.  5.50  to  5.52  show  the  results  of  the 
comparisons. 

FNOC  climatology  was  in  reasonable  agreement  with  the 
1982  mean  profile  (assumed  non-El  Nino  year),  Fig  5.50  .  The 
major  discrepancy  between  climatological  and  measured 
profiles  occurred  in  the  upper  part  of  the  water  column 
(surface  to  50m)   which  could  be   expected  since  this  is  the 
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region  subject  to  maximum  atmospheric  forcing  and,  there- 
fore, the  short  time  scale  variability  impossible  to  predict 
with  climatology. 

The  presence  of  a  thermal  anomaly  in  the  depth/ 
temperature  profile  was  clearly  evident  when  the  1982 
profile  was  compared  to  that  from  1983,  Fig.  5.51.  The 
assumed  El  Nino  event  was  present  as  a  layer  of  warmer 
waters  between  the  near  surface  and  approximately  200m  with 
a  maximum  anomaly  of  2. 5°C  at  a  depth  of  60m.  The  inability 
of  a  static  climatology  to  represent  such  an  event  was 
demonstrated  by  comparing  the  1983  mean  depth/temperature 
profile  with  that  obtained  from  FNOC  climatology,  Fig.  5.52. 

To  obtain  a  measure  of  the  operational  significance  of 
an  anomalous  event  such  as  the  one  described  above  transmis- 
sion loss  curves  were  computed  from  the  FACT9H  transmission 
loss  model  for  theJuly  1983  climatological  and  mean  depth/ 
temperature  profiles  .  The  only  measurable  effect,  Figs. 
5.53  to  5.55  was  a  reduction  in  convergence  zone  range  in 
the  curves  derived  ficm  the  climatological  profiles  (approx- 
imately 3%)  and  a  decrease  in  convergence  zone  gain 
(approximately  5db)  . 

C.   RESULTS  OF  RECENT  FNOC  STUDIES 

Recent  modifications  have  been  made  to  both  EOTS  and 
TEOTS  analysis  models  (Frost,  pers.  comm.)  and,  in  partic- 
ular, the  model  sensitivity  to  input  changes  in  P1D  (MLD) 
has  been  modified.  The  incorporation  of  the  TEOTS  analysis 
into  the  finer  grid  regional  model  has  also  been  achieved 
for  certain  areas.  The  effect  of  these  modifications  is  the 
subject  of  an  ongoing  investigation  by  FNOC,  based  on  inputs 
from  the  Bay  of  Biscay  and  the  Western  Mediterranean  Sea 
regional  analyses. 
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The  initial  results  from  the  study  are  promising  with 
TEOTS  apparently  showing  considerable  skill  in  predicting 
changes  in  layer  depth  even  in  the  absence  of  input  data 
(Frost,  pers.  com.)  .  No  information  is  available  from  this 
study  related  to  the  skill  of  the  analysis  in  resolving  the 
horizontal  structure  in  the  region. 
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Figure   5.1        Contoured   SST   Fields   for   EOTS/Cruise   AI 

(contour    units  °Cx10) . 
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b.  EOTS  Analysis  16  Jun 
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Figure  5.2 


Contoured  HID  Fields  for  EOTS/Cruise  AI 
(contour  units  m)  . 
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Figure   5,3        Contoured   BLG   Fields  for   EOTS/Cruise   AI 
(contour   units  oc/1000m) . 
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Figure   5.4        Contoured   THG   Fields  for   EOTS/Cruise   AI 
(contour   units   oc/1000m). 
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Figure   5.5        Contoured  SST   Fields   for  TEOTS/Cruise  AI 

(contour   units   °Cx10) . 
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Figure  5.6 


Contoured  MLD  Fields  for  TEOTS/Cruise  AI 
(contour  units  a)  . 
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Figure  5.7    Contoured  BLG  Fields  for  TEOTS/Cruise  AI 
(contour  units  oc/1000m) . 
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Figure   5-8 


Contoured  THG   Fields   for  TEOTS/Cruise   AI 
(contour  units  <>C/1000m)  . 
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Figure  5.9        Contoured   SST   Fields   for  EOTS/Cruise   All 

(contour    units   °Cx10) . 
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a.  Cruise  All  Analysis 
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Figure  5.10        Contoured   HLD   Fields   for  EOTS/Cruise   All 

(contour   units   m)  . 
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Figure   5.11        Contoured   BLG   Fields  for   EOTS/Cruise   All 
(contour   units  oc/1000m) . 
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Figure   5.12        Contoured  THG   Fields  for   EOTS/Cruise   All 


(contour   units   oc/1000m) 
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Figure   5.13        Contoured   SST   Fields   for   TEOTS/Cruise   All 

(contour    units   °Cx10) . 
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Figure  5. 14 


Contoured  HLD  Fields  for  TEOTS/Cruise  All 
(contour  units  a) . 
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Figure  5-15        Contoured   BLG   Fields  for  TEOTS/Cruise   All 

(contour  units   oc/1000) . 
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Figure  5.16 


Contoured  THG  Fields  for  TEOTS/Cruise  All 
(contour  units  oc/1000)  . 
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Figure  5.  17 


Contoured  SST  Fields  for  EOTS/Cruise  AIII 
(contour  units  °Cx10) . 
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Figure   5.18        Contoured   BLD   Fields   for   EOTS/Cruise   AIII 

(contour   units   n)  . 
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Figure  5-19        Contoured   BLG   Fields   for   EOTS/Cruise   AIII 

(contour   units    °C/1000) . 
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Figure  5.20        Contoured  TBG  Fields   for  EOTS/Cruise  AIII 

(contour   units   «C/1000) . 
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Figure  5.21        Contoured   SST   Fields   for   TEOTS/Cruise   AIII 

(contour   units   °Cx10) . 
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Figure   5.22        Contoured   HLD   Fields  for  TEOTS/Cruise   AIII 

(contour   units   m)  . 
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Figure  5.23        Contoured   BLG   Fields   for  TEOTS/Cruise   AIII 

(contour   units   <>C/1000)  . 
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Figure   5.24        Contoured  THG  Pields   for  TEOTS/Cruise   AIII 

(contour   units   oc/1000) . 
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Figure  5.25        SST    Difference  Distribution   EOTS/Cruise  AI 
a.    16   Jun   Analysis   b.    23   Jun   Analysis. 
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Figure   5.26        HLD    Difference   Distribution   EOTS/Cruise  AI 
a.    16   Jun   Analysis   b.    23   Jun   Analysis. 
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Figure  5.27        BLG    Difference   Distribution   EOTS/Cruise  AI 
a.    16   Jun   Analysis   b.    23   Jun   Analysis. 
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Figure  5.28        THG    Difference  Distribution   EOTS/Cruise  AI 
a.    16   Jun   Analysis   b.    23   Jun   Analysis. 
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Figure  5.29    SST  Difference  Distribution  TEOTS/Cruise  AI 
a.  16  Jun  Analysis  b.  23  Jun  Analysis. 
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Figure   5.30        HLD   Difference   Distribution   TEOTS/Cruise  AI 
a.    16   Jun   Analysis   b.    23   Jun   Analysis. 
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Figure  5.31    BLG  Difference  Distribution  TEOTS/Cruise  AI 
a.  16  Jun  Analysis  b.  23  Jun  Analysis. 
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Figure  5.32   THG  Difference  Distribution  TEOTS/Cruise  AI 
a.  16  Jun  Analysis  b.  23  Jun  Analysis. 
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Figure   5-33        SST   Difference   Distribution   EOTS/Cruise   All 
a.    30   Jun   Analysis   b.    7  Jul   Analysis. 
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Figure   5.34        HLD   Difference   Distribution    EOTS/Cruise   All 
a.    30   Jun   Analysis   b.    7  Jul   Analysis. 
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Figure   5-35        BLG   Difference   Distribution   EOTS/Cruise   All 
a.    30   Jun   Analysis  b.    7   Jul   Analysis. 
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Figure    5.36        THG   Difference   Distribution   EOTS/Cruise   All 
a.    30   Jun   Analysis  b.    7   Jul   Analysis. 
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Figure  5.37   SST  Difference  Distribution  TEOTS/Cruise  All 
a.  30  Jun  Analysis  b.  7  Jul  Analysis. 
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Figure  5.38   MLD  Difference  Distribution  TEOTS/Cruise  All 
a.  30  Jun  Analysis  b.  7  Jul  Analysis. 
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Figure  5-39   BLG  Difference  Distribution  TEOTS/Cruise  All 
a.  30  Jun  Analysis  b.  7  Jul  Analysis. 
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Figure   5.40        THG   Difference   Distribution   TEOTS/Cruise  All 
a.    30   Jun   Analysis   b.    7   Jul   Analysis. 
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Figure  5.41 


SST  Difference  Distribution  EOTS/Cruise  All 
10  Jul  Analysis. 
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Figure  5.42 


MLD  Difference  Distribution  EOTS/Cruise  III 
10  Jul  Analysis. 
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Figure  5.43 


BLG  Difference  Distribution  EOTS/Cruise  All 
10  Jul  Analysis. 
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Figure  5.44 


THG  Difference  Distribution  EOTS/Cruise  All 
10  Jul  Analysis. 
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Figure   5-45 


SST  Difference  Distribution  TEOTS/Cruise  All 
10  Jul  Analysis. 
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Figure  5.46 


MLD  Difference  Distribution  TEOTS/Cruise  All 
10  Jul  Analysis. 
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Figure   5.47 


BLG  Difference  Distribution  TEOTS/Cruise  All 
10  Jul  Analysis. 
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Figure  5.48 


THG  Difference  Distribution  TEOTS/Cruise  All 
10  Jul  Analysis. 
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Figure  5. 49    SST  Contour  Fields:  (a)  EOTS  23  Jun  Analysis 
(b)  Author* s  Contouring  Program 
(same  inputs  as  EOTS  23  Jun  Analysis) . 
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Figure   5.50        Temperature   Anomaly   vs.    Depth 
Jul    1982-FNOC  Climatology. 
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Figure   5.51 


Temperature  Anomaly  vs.  Depth 
Jul  1983-Jul  1982. 
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Figure   5.52        Temperature  Anomaly  vs.    Depth 
Jul    1983-FHOC  Climatology. 
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Figure  5-53        Transmission   loss  Curves-S/S   Disposition 
a.    Climatology  b.    0PT0MA5    AIII    (mean   profile). 
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Figure  5.54    Transmission  loss  Curves-S/D  Disposition 
a.  Climatology  b.  0PT0MA5  AIII  (mean  profile) . 
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Figure  5.55    Transmission  Loss  Curves-D/D  Disposition 
a.  Climatology  b.  0PT0MA5  AIII  (mean  profile) . 
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PARAMETER 


SST 
(°  C) 


TABLE    II 
Comparison   Statistics   Cruise   AI 


MODEL 


ANALYSIS  DATE 


mean 


std.    deviation 


median 


skevmess 


EOTS 


16   Jun 


-0.80 


0.50 


-0.60 


-0.40 


23   Jun 


0.10 


0.30 


0.20 


-0.30 


TEOTS 


16   Jun 


-0.90 


0.50 


-1.00 


-0.30 


23   Jun 


-0.50 
0.40 


-0.60 


-0.20 


MLD 
(m) 


mean 


-19.4 


-17.9 


-6.50 


std.    deviation 


7.70 


7.30 


.70 


median 


-15.0 


•15.0 


-5.00 


skevmess 


0.10 


0.20 


-0.10 


-5.70 


7.70 


-5.00 


-0.20 


BLG 
(°C/100m) 


mean 


-1.26 


-0.05 


-3.24 


std.    deviation 


1.16 


0.70 


0.90 


med  ian 


-0.50 


-0.50 


-3.50 


skevmess 


-0.80 


0.50 


-0.10 


-1.60 


1.00 


-1.50 


-0.40 


THG 
(°C/100ra) 


mean 


On  r\ 
.'JO 


0.02 


-0.11 


std.    deviation 


median 


skevmess 


o.i: 


0.05 


0.50 


0.07 


0.05 


0.30 


0.12 


-0.05 


-0.60 


-0.09 


0.09 


-0.05 


-0.60 
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TABLE  III 
Comparison  Statistics  Cruise  All 

PARAMETER 

MODEL 

EOTS 

TEOTS 

ANALYSIS  DATE 

30  Jun 

07  Jul 

30  Jun 

07  Jul 

SST 
(°C) 

mean 

0.90 

1.03 

0.40 

-0.04 

std.  deviation 

0.40 

0.50 

0.60 

0,60 

median 

1.00 

1.00 

0.60 

0.20 

skevmess 

-0.90 

-0.70 

-0.02 

0.01 

MLD 
(m) 

mean 

-30.9 

-26.6 

-11,2 

-17.4 

std.  deviation 

8,30 

7.70 

16.8 

7,89 

median 

-25.0 

-25.0 

-15.0 

-15.0 

skewness 

-0.20 

-0.60 

2.16 

-0.80 

BLG 
("C/lOOm) 

me.  in 

-1.20 

-0.30 

-2.80 

-4.24 

std.  deviation 

1,90 

1.40 

1.60 

1,49 

median 

-1.50 

-0.50 

-3.50 

-4.50 

skewness 

0.90 

1.40 

1.20 

0.90 

THG 
fC/lOOm) 

mean 

0.09 

0.25 

0.08 

-0.04 

std.  deviation 

0.07 

0.12 

0.11 

0.11 

median 

0.20 

0.30 

0.15 

-0,05 

skewness 

-0.60 

-1.00 

-0.30 

-0.20 
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TABLE   17 
Comparison   Statistics  Cruise 

— , .              1 

AIII 

PARAMETER 

MODEL 

EOTS 

TEOTS 

ANALYSIS   DATE 

10  Jul 

21 

Jul 

10   Jul 

21   Jul 

SST 

(°  c) 

mean 

0.50 

0.12 

std.    deviation 

0.60 

0.60 

median 

0.60 

0.20 

skewness 

0.20 

0.10 

MLD 
(m) 

mean 

-5.60 

1.50 

std.    deviation 

16.90 

13.10 

median 

-5.00 

5.00 

skewness 

-0.40 

-0.70 

BLG 
(*C/100m) 

mean 

0.05 

-3.90 

std.    deviation 

1.33 

1.15 

median 

0.50 

-3.50 

skewness 

0.20 

-0.05 

THG 
(°C/100m) 

mean 

0.20 

0.04 

std.    deviation 

0.14 

0.14 

median 

0.20 

0.05 

skewness 

-0.20 

0.30 



124 


VI.     ACOUSTIC    I HP ACT    OF   EBBOBS    IN    FNOC    OCEAN    THEBMAL 

STBUCTUBE    ANALYSIS 

A.  DEFINITION  OF  ACOUSTIC  TEBHS 

•  Convergence  Zone  (CZ)  -  Region  of  focussing  or  convergence 
cf   sound. 

•  Direct  Range  -  The  shortest  distance  from  the  receiver  at 
which    the   transmission   loss    exceeds    the    Figure   of    Merit. 

•  Figure  of  Merit  (POM)  -  The  maximum  allowable  one  way 
transmission  loss  (passive  sonars)  for  which  a  50%  prob- 
ability  of   detection    is  possible    for    a   particular   sonar. 

•  Leakage  -  Loss  of  energy  from  a  surface  duct  to  the  main 
sound  channel. 

•  Low-Frequency  Cut-Off  -  Frequency  below  which  no  trapped 
normal   mode   may   exist   for   a    given  channel    or   duct. 

•  Sound-Speed  Profile  (SSP)  -  Depth/sound-speed  combinations 
defining  the  vertical  variation   in   sound   speed. 

•  Transmission  Loss  (TL)  -  Acoustic  parameter  which  quanti- 
tatively describes  the  weakening  of  sound  intensity  between 
a  point  1m  from  the  source  and  a  point  at  a  distance  in  the 
ocean. 

•  Vertical  Beam-Width  -  Included  angle  between  the  maximum 
upward   and   downward    propagating  rays. 

B.  ASSUMPTIONS 

The  cnly  variable  in  the  range  calculations  was  to  be 
the  SSP.  To  achieve  this  goal  certain  representative, 
acoustic  parameters  were  assumed  constant  throughout  and 
were  used  as  standard  inputs  to  all  transmission  loss 
calculations. 

•  The  ocean  bottom  was  assumed  flat. 
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•  A  standard  bottom- less/grazing  angle  curve  applied  to  all 
frequencies. 

•  A    flat   sea  surface    (sea  state    1)       input   to   the    FACT   models 
(assumed   flat   in   the    FE   model) 

•  The  FOM's  used  in  the  range  calculations  were  80db  and 
90db.  These  figures  are  representative  of  the  range  of  FOM»s 
encountered   in   operational  systems. 

•  Source/receiver  dispositions   used    were: 

20m/20m-shallow/shallow     (S/S) 

20m/200m-shallow/deep    (S/D) 

200m/200m-deep/deep  (D/D) 
The  S/S  disposition  was  intended  to  give  a  measure  of 
the  surface  duct  acoustic  propagation;  S/D  of  cross-layer 
propagation  and  D/D  of  deep  sound  channel  propagation.  All 
three  combinations  can  be  utilized  by  presently  operational, 
shipberne   systems. 

C.       DERIVATION    OF   SOUBD-SPEED    PROFILES 

Depth/temperature  profiles  were  constructed  for  both  the 
EOTS  and  0PT0MA5  analyses  at  the  positions  of  maximum  error, 
between  the  two  analyses,  for  SST,  MLD,  BLG  and  THG.  The 
position  of  the  maximum  error  was  obtained  from  the  differ- 
ence field  discussed  in  Chapter  V  and  the  depth/temperature 
profiles  were  generated  by  utilizing  the  contour  fields 
derived    as   described    in  Chapter  IV. 

SSP's  were  computed  from  the  temperature  profiles  using 
Wilson's  Equation.  The  depth/salinity  information  was 
obtained   from   the    DE   STEIGEDR,    deep    CTD    data    (Chapter   II) 

Only  the  data  relevant  to  cruises  AI  and  All  were 
utilized.  Long  range  prediction  (CZ)  comparisons  were  made 
between  profiles  based  on  the  0PT0MA5,  EOTS  and  TEOTS  anal- 
yses. Direct  range  comparisons  were  limited  to  profiles 
based   on    the  0PT0MA5    and   EOTS    analyses   only. 
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The  SSP's  derived  for  each  model  analysis  (and  all 
subsequent  acoustic  calculations  based  on  those  SSP's)  were 
designated  according  to  the  parameter,  and  analysis  date  eg. 
SST1  refers  to  a  profiles  derived  at  the  position  of  maximum 
difference  between  modeled  and  0PT0MA5  SST  for  the  initial 
model  analysis  of  the  particular  0PT0MA5  cruise  being 
studied;  MLD2  refers  to  the  profile  derived  at  the  position 
of  maximum  difference  in  MLD's  for  the  final  model  analysis 
of  the  cruise,  etc. 

D.   TBANSMISSION  LOSS  MODEL  COMPARISONS 

Differences  in  model  physics  made  direct  comparisons 
between  FACT9H  and  the  PE  Models  unwieldy,  as  did  limita- 
tions in  the  model  computer  programs.  One  of  the  major 
differences  was  that  the  PE  model  has  a  vertical  beamwidth 
limitation  which  causes  a  program  termination  if  exceeded. 
The  limitation  is  a  function  of  the  SSP  and,  for  the 
profiles  utilized  in  this  study,  the  maximum  vertical  beam- 
width  which  could  be  maintained  was  20°.  No  such  limitation 
exists  fcr  the  FACT9H  model  which  sums  acoustic  intensities 
from  all  forward  directions,  equivalent  to  a  vertical  team- 
width  of  180°.  The  effect  of  this  beamwidth  limitation  is 
that  the  PE  model  will  predict  consistantly  greater  trans- 
mission loss  for  a  given  range  than  will  the  FACT9H  model, 
assuming  all  other  inputs  to  each  model  are  identical. 
However,  the  effect  of  SSP  changes  relating  to  temperature 
profile  differences  between  modeled  and  OPTOMA5  ocean 
thermal  structures  were  computed  on  a  model-by-model  basis 
and  the  results  compared,  Table  V.  An  example  of  the  trans- 
mission loss  curves  computed  from  the  two  models  for  the 
same  SSI,  Fig.  6.1,  clearly  shows  the  effect  of  the 
differing  model  physics. 
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Table  V  shows  the  Parameter  Error  (difference  between 
OPT0MA5  and  EOTS  values  for  the  parameter,  expressed  as  a 
percentage  of  the  0PTCMA5  value) ;  the  Range (E)  (direct  range 
(km)  computed  from  the  transmission  loss  curve  based  on  the 
0PT0MA5  data) ;  Range  Error  (E)  (difference  between  direct 
range  obtained  utili2ing  0PT0MA5  data  and  that  utilizing 
EOTS  analysis  data,  expressed  as  a  percentage  of  R).  A 
negative  sign  indicates  that  the  model  analysis  had  produced 
an  overestimate  of  that  particular  parameter.  The  compari- 
sons shown  in  Table  V  are  for  a  600Hz  frequency. 

The  lack  of  agreement  between  the  two  transmission  loss 
models  was  apparent  with  the  PE  model  predicting  signifi- 
cantly shorter  ranges  than  the  FACT9H  model.  The  PE  model 
was  also  far  more  sensitive  to  variations  in  the  input  SSP, 
particularly  when  both  source  and  receiver  were  shallow,  and 
showed  a  dependence  on  the  SSP  for  the  S/D  and  D/D  disposi- 
tions, which,  as  discussed  in  later  sections,  was  not  the 
case  for  the  FACT9H  model.  Convergence  zone  ranges  computed 
by  both  FACT9H  and  the  PE  model  were  in  good  agreement. 
However,  the  gain  in  acoustic  energy  at  the  convergence 
zones  was  approximately  5db  lower  in  the  PE  model  computa- 
tions than  those  of  FACT9H. 

The  remaining  comparisons  discussed  in  this  chapter  were 
based  on  the  FACT9H  transmission  loss  computations. 

E.   SUHHARY  OF  ACOUSTIC  COMPUTATIONS 

When  comparing  ranges  predicted  utilizing  measured  and 
modeled  profiles,  an  error  of  less  than  20%  was  assumed 
operationally  insignificant.  This  assumption  was  telieved 
reasonable  when  the  uncertainties  associated  with  the  FACT9H 
acoustic  model  accuracy  (eg. ,  the  lack  of  attenuation  due  to 
absorption  and  volumetric  scattering)  are  considered. 
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1 .   Lon£  Range  (CZ)  Predictions 

Comparisons  of  CZ  ranges  computed  using  SSP's 
derived  from  OPTOMA5,  EOTS  and  TEOTS  data,  Tables  VI  to  IX 
revealed  a  relative  lack  of  sensitivity  of  the  CZ  range  to 
thermal  structure  variations  of  the  order  of  magnitude 
present  between  modeled  and  sea  truth  analyses.  The  average 
CZ  range  computed  from  sea  truth  data  was  50km,  varying 
between  44km  and  62km.  Ranges  computed  utilizing  SSP's 
derived  from  FNOC  model  analyses  were  within  10??  cf  those 
computed  utilizing  sea  truth  data.  The  maximum  errors  were 
observed  when  both  source  and  receiver  were  shallow. 

2  •   Direct  R anc[e  Compear isons 

Subjective  analysis  of  the  transmission  loss  curves 
computed  for  the  long  range  comparisons  discussed  above, 
combined  with  an  understanding  of  ray-theory  physics,  led  to 
the  direct  range  comparisons  being  concentrated  on  the 
surface  duct  propagation  paths  (i.e.,  shallow  source  and 
receiver).  A  limited  number  of  comparisons  were  made  between 
direct  ranges  obtained  from  the  transmission  loss  curves 
computed  utilizing  OPTOMA5  and  those  utilizing  EOTS  derived 
profiles  for  cross-layer  and  below  layer  propagation  paths. 
The  results  from  these  comparisons,  Tables  X  and  XI  confirm 
the  supposition  that  the  FACT9H  computed  transmission  losses 
for  the  S/D  and  D/D  propagation  paths  are  relatively  insen- 
sitive to  SSP  variations  of  the  order  of  magnitude  present 
between  profiles  derived  from  0PT0MA5  data  and  those  from 
the  EOTS  analyses.  In  the  remainder  of  this  discussion 
'range  error1  will  refer  to  the  difference  between  ranges 
calculated  from  transmission  loss  curves  computed  from  sea 
truth  profiles  and  those  computed  from  EOTS  profiles. 

The  maximum  range  error  for  the  cross- layer  propaga- 
tion path  was  18%   on  a  range  of  8km,   with   the  majority  of 
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the  range  errors  within  10%  of  a  sea  truth  range  of  8km  at 
the  50Hz  frequency.  However,  a  weak  frequency  dependence  was 
observed  and  errors  cf  the  order  of  15%  on  a  range  of  8km 
were  measured  at  600Hz,  associated  with  the  lower  FOM.  The 
maximum  error  also  occurred  at  this  frequency,  associated 
with  the  lower  FOM,  Table  X  - 

The  range  errors  measured  for  the  below  layer  propa- 
gation paths  were  generally  negligible  (<5%)  and  independent 
of  frequency  and  FOM,  Table  XI  .  However,  for  profiles 
SSTIand  SST2  errors  of  the  order  of  15%  on  a  range  of  18km 
for  the  lower  FOM  and  15%  on  a  range  of  77km  for  the  higher 
FOM,  were  measured.  The  presence  of  a  weak  shallow  sound 
channel  in  the  sea  truth  SSP  is  believed  to  be  the  cause  of 
these  significantly  higher  range  errors  and  would  also 
explain  the  longer  direct  ranges  associated  with  the  sea 
truth  profile  (i.e. ,  18km  compared  to  the  average  of  8km  for 
the  remaining  profiles) 

For  the  surface  duct  propagation  paths  (S/S) ,  trans- 
mission loss  curves  were  computed  at  frequencies  of  50,  300 
600  and  1000Hz  for  profiles  associated  with  cruise  AI  and  at 
300,  600  and  1000Hz  for  those  associated  with  cruise  All. 
The  50Hz  frequency  was  removed  from  the  comparisons  in 
cruise  All,  Table  XIII,  once  the  range  errors  were  shown  to 
be  negligible  at  this  frequency  for  the  cruise  AI  compari- 
sons, Table  XII  . 

From  the  range  comparisons.  Table  XII  and  XIII  it 
can  be  seen  that  at  the  higher  FOM  and  for  frequencies  below 
1000Hz,  the  transmission  loss  is  insensitive  to  SSP  varia- 
tions of  the  order  of  magnitude  observed  with  a  range  error 
of  less  than  10%  on  a  mean  range  of  55km. 

For  range  calculations  based  on  the  lower  FOM, 
significant  range  errors  were  observed  for  frequencies  of 
300Hz  and  above.  The  range  error  was  also  observed  to  be  a 
function  of  frequency.  Plots  of  range  error  versus  parameter 
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error  for  the  four  parameters  being  studied  revealed  a 
strong  dependence  of  the  range  error  on  the  MLD  error,  Fig. 
6.2  .  The  curves  shown  in  Fig.  6.2  were  fitted  to  the  points 
in  an  RMS  manner  and  do  not  imply  a  functional  relationship 
between  MID  and  range  error  but  are  intended  to  shew  a 
general  trend.  The  scatter  of  the  points  in  Fig.  6.2  indi- 
cate the  complex  interaction  between  the  various  parameters. 
From  the  curves  it  can  be  seen  that  the  range  error 
increased  with  increasing  MLD  error  and  frequency,  for 
freguecies  above  approximately  300Hz.  The  presence  cf  a 
non-zero  range  error  for  a  zero  MLD  error  was  a  further 
indication  of  the  effect  of  the  other  parameters.  Plcts  of 
range  error  versus  errors  in  the  remaining  three  parameters 
were  also  generated,  Figs.  6.3  to  6.5  These  curves  appeared 
to  indicate  a  marked  dependence  of  range  error  on  each 
parameter.  However,  tlots  of  MLD  error  versus  the  errors  in 
the  remaining  parameters,  Figs.  6.6  to  6.8  revealed  that  the 
major  range  error  dependence  was  indeed  on  MLD  error;  e.g. 
the  trend  in  range  error  when  compared  to  SST  error  is  in 
the  same  general  sense  as  the  trend  in  MLD  error  when  that 
too  is  compared  to  SSI  error. 
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Figure  6.3        Plot  of   SST   Error   vs.    Range  Error 
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Figure  6.5   Plot  of  THG  Error  vs.  Range  Error 
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Figure  6.6    Plot  of  SST  Error  vs.  MID  Error. 
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Figure  6.7    Plot  of  BIG  Error  vs.  HLD  Error, 


138 


u 
a: 

Q. 


2 

<3 


1               1 
-5 

1 

i 

5 

10 

• 

15 

• 

• 

• 

o 
o 

CM 

1 

i 

o 
o 

1 

. 

o 
o 

m 

• 

1 

ATHG  (PERCENT) 


Figure  6.8    Plot  of  THG  Error  vs.  MLD  Error, 
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VII.  SDHHARY  AND  CO NC LOS IONS 

The  grid  spacing  used  in  the  FNOC  ocean  thermal  analysis 
models  (i.e.,  40km  for  EOTS  and  320km  for  TEOTS)  made  it 
extremely  unlikely  that  they  would  be  capable  of  recovering 
the  mesoscale,  horizontal  variability  (of  the  order  of  30  to 
40km)  dominating  the  study  domain. 

Within  the  limitations  of  the  grid  spacing  used,  the 
EOTS  SST  analysis  responded  well  to  input  data.  However,  the 
analysis  methodology  is  believed  to  have  smoothed  the  input 
field  to  the  extent  that  significant  detail  was  lost. 

The  hemispherical  TEOTS  SST  analysis  responded  to  data 
input  in  magnitude  but  was  totally  insensitive  to  the  hori- 
zontal shape  information  provided  by  the  input  data.  This 
lack  of  response  is  believed  to  be  a  function  of  the  coarse 
grid  spacing  of  the  TEOTS  analysis  field. 

Both  EOTS  and  TEOTS  were  insensitive  to  input  mixed 
layer  depth  (MLD)  information.  The  EOTS  MLD  analysis 
remained  almost  constant  throughout  the  study  period  in  both 
magnitude  and  shape.  The  TEOTS  MLD  analyses  did  change  over 
the  study  period  and  were,  on  average,  in  closer  agreement 
with  the  measured  MLD  magnitudes,  but  had  an  almost  constant 
shape  field.  The  lack  of  detail  in  both  EOTS  and  TEOTS  MLD 
analyses  was  in  marked  contrast  to  the  complex  and  variable 
measured  MLD  shape  fields.  The  improved  performance  of  the 
TEOTS  analysis  over  the  EOTS  analysis  is  believed  to  be  the 
result  of  the  influence  of  the  TOPS  element  on  the  analysis. 

The  response  of  the  models  to  data  input  for  depths 
below  the  MLD  was  inconclusive  since  the  parameters  tested 
(BLG  and  THG)  were  functions  of  the  MLD  and,  thus,  errors  in 
the  MLD  analysis  would  propagate  into  these  remaining 
fields. 
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The  FNOC  procedure  of  using  a  single  climatological 
depth/temperature  profile  for  areas  such  as  the  study  domain 
was  shown  to  be  inaccurate.  The  study  domain  was  shown  to 
vary  on  a  secular  time  scale  (El  Nino  effect)  with  a  2  to 
3°C  thermal  anomaly,  compared  to  FNOC  climatology,  measured 
at  depths  of  50  to  100m  during  the  study  period.  The  effect 
of  a  poor  climatology  would  be  reflected  in  the  model  anal- 
yses in  the  absence  of  input  data. 

For  acoustic  frequencies  below  1000Hz  and  based  on 
results  from  the  FACT9H  Transmission  Loss  model,  the  errors 
in  the  EOTS  and  TE01S  analyses  were  insignificant  for  (a) 
long  range  (convergence  zone)  predictions,  (b)  direct  range 
predictions  for  sonar  systems  with  a  figure  of  merit  of  90db 
or  greater,  and  (c)  direct  range  predictions  for  cross  layer 
and  below  layer  source/receiver  dispositions.  The  errors  in 
the  model  analyses  were  found  to  be  significant  for  sonar 
systems  with  a  figure  of  merit  of  80db  or  less  when  both 
source  and  receiver  were  shallow. 

The  errors  in  modeled  mixed  layer  depth  were  found  to 
have  the  greatest  effect  on  acoustic  prediction  errors,  and 
the  effect  was  found  to  be  directly  proportional  to  acoustic 
frequency. 

A  comparison  between  the  FACT9H  and  PE  Transmission  Loss 
models  revealed  significant  differences  between  the  outputs 
of  the  two  models. 

The  PE  model  appeared  far  more  sensitive  to  errors  in 
the  temperature  profile  than  did  the  FACT9H  model.  The 
direct  range  errors  obtained  from  the  PE  model  outputs, 
based  on  modeled  and  measured  temperature  profiles,  showed 
no  dependence  on  the  figure  of  merit  and  less  dependence  on 
the  source/receiver  disposition  than  did  those  based  on  the 
FACT9B  model  outputs. 

The  disparity  between  FACT9H  and  PE  model  sensitivities 
made   an   absolute  determination   of   the   operational 
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significance  of  errors  in  the  EOTS  and  TEOTS  analyses 
impossible.  However,  the  maximum  effect  of  depth/temperature 
profile  differences  on  direct  range  predictions  was 
significant,  whichever  transmission  loss  model  was  used, 
when  both  source  and  receiver  were  shallow. 
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VIII-  RECOMMENDATIONS 

The  staff  at  FNOC  are  presently  restructuring  the  ocean 
thermal  analysis  model  algorithms  to  enable  the  models  to 
respond  more  readily  to  input  data,  particularly  the  PLD 
segment  cf  the  analysis.  The  physics  of  the  TEOTS  analysis 
are  also  being  incorporated  into  the  finer  grid,  regional 
analyses.  The  effectiveness  of  such  changes  should  be  tested 
in  a  similar  areal  manner  to  that  of  the  present  study,  and 
for  many  more  areas  of  the  world's  ocean,  to  establish  the 
operational  value  of  such  an  analysis  system.  The  minimum 
data  input  required  to  give  a  representative  analysis  under 
varying  conditions  should  also  be  established,  in  order  that 
assets  may  be  deployed  effectively  and  economically. 

For  regions  where  anomalous  events  are  known  to  occur, 
the  feasability  of  utilizing  a  quasi-dynamic  climatology, 
based  on  hindcasting,  should  be  studied,  i.e.,  climatology 
could  he  updated  when  a  persistent  anomalous  event  is  known 
to  be  present.  This  would  give  a  more  representative  anal- 
ysis, at  depths  below  those  influenced  by  the  TOPS  elements 
of  the  TECTS  analysis,  in  the  absence  of  input  data. 

Satellite  IE  imagery  could  be  a  powerful  asset  in  data- 
sparse  areas  of  the  oceans,  providing  a  necessary  input  to 
ocean  thermal  analysis  models.  However,  the  results  of  the 
present  study  have  shown  that  the  relationship  between  SST 
and  sub-surface  features  is  not  always  an  obvious  one. 
Studies  should  be  undertaken  to  relate  surface  features  with 
sub-surface  structure,  either  through  direct  correlations  or 
modeling.  Such  relationships  would  also  assist  the  planning 
of  sampling  strategies  aimed  at  achieving  the  most  econom- 
ical and  efficient  use  of  assets  whilst  providing  sufficient 
input  to  the  ocean  thermal  analysis  model  to  achieve  the 
desired  accuracy. 
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The  effect  on  the  higher  frequency ; i. e- ,  greater  than 
1Khz,  active  sonar  performance  of  errors  in  the  model  anal- 
ysis fields  should  be  examined.  The  loss  of  vertical  struc- 
ture could  prove  important.  In  particular,  areas  where 
significant  vertical  features  are  present;  e.g.,  shallow 
sound  channels,  should  be  examined  and  the  feasibility  of  an 
additional  floating  level  in  the  analysis,  to  represent  the 
depth  of  such  features,  should  be  addressed. 

The  effect  of  mesoscale  variations,  on  scales  such  as 
those  observed  in  the  present  study,  on  the  other  parameters 
in  the  sonar  equation  should  be  examined;  e.g.,  ambient 
noise  changes  due  to  changes  in  the  marine  ecosystem  in  the 
region  of  the  weak  ocean  fronts  associated  with  upwelling 
regions. 

.  Operationally  important  areas  are,  in  general,  the  more 
densely  sampled  regions  of  the  ocean.  However,  the  present 
fixed  regional  boundaries  utilized  in  the  FNOC  analysis 
models  are  only  broadly  representative  of  these  operational 
areas  and  generally  extend  beyond  the  limits  of  prime 
interest.  As  a  consequence,  the  shape  field  used  in  the 
analysis  is  influenced  by  possibly  less  well  sampled  areas 
outside  the  region  of  interest  with  possible  erroneous 
smoothing  effects  on  the  horizontal  and  vertical  thermal 
structure  that  the  analysis  would  otherwise  produce.  To 
remove  such  a  possibility,  the  utilization  of  a  flexible 
regional  boundary,  fitted  to  the  area  of  interest  ,  and  with 
a  grid  spacing  reprasentative  of  the  horizontal  variability 
of  the  region,  should  be  studied.  Such  an  approach  as  the 
one  outlined  above  would  also  facilitate  the  incorporation 
of  a  more  detailed,  possibly  dynamic,  climatology  into  the 
analysis. 
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